Portability of Scientific Workflows in NGS Data Analysis: A Case Study

The analysis of next-generation sequencing (NGS) data requires complex computational workflows consisting of dozens of autonomously developed yet interdependent processing steps. Whenever large amounts of data need to be processed, these workflows must be executed on a parallel and/or distributed systems to ensure reasonable runtime. Porting a workflow developed for a particular system on a particular hardware infrastructure to another system or to another infrastructure is non-trivial, which poses a major impediment to the scientific necessities of workflow reproducibility and workflow reusability. In this work, we describe our efforts to port a state-of-the-art workflow for the detection of specific variants in whole-exome sequencing of mice. The workflow originally was developed in the scientific workflow system snakemake for execution on a high-performance cluster controlled by Sun Grid Engine. In the project, we ported it to the scientific workflow system SaasFee that can execute workflows on (multi-core) stand-alone servers or on clusters of arbitrary sizes using the Hadoop. The purpose of this port was that also owners of low-cost hardware infrastructures, for which Hadoop was made for, become able to use the workflow. Although both the source and the target system are called scientific workflow systems, they differ in numerous aspects, ranging from the workflow languages to the scheduling mechanisms and the file access interfaces. These differences resulted in various problems, some expected and more unexpected, that had to be resolved before the workflow could be run with equal semantics. As a side-effect, we also report cost/runtime ratios for a state-of-the-art NGS workflow on very different hardware platforms: A comparably cheap stand-alone server (80 threads), a mid-cost, mid-sized cluster (552 threads), and a high-end HPC system (3784 threads)

Calibration of TCCON column-averaged CO2: the first aircraft campaign over European TCCON sites

The Total Carbon Column Observing Network (TCCON) is a ground-based network of Fourier Transform Spectrometer (FTS) sites around the globe, where the column abundances of CO2, CH4, N2O, CO and O2 are measured. CO2 is constrained with a precision better than 0.25% (1-σ). To achieve a similarly high accuracy, calibration to World Meteorological Organization (WMO) standards is required. This paper introduces the first aircraft calibration campaign of five European TCCON sites and a mobile FTS instrument. A series of WMO standards in-situ profiles were obtained over European TCCON sites via aircraft and compared with retrievals of CO2 column amounts from the TCCON instruments. The results of the campaign show that the FTS measurements are consistently biased 1.1% ± 0.2% low with respect to WMO standards, in agreement with previous TCCON calibration campaigns. The standard a priori profile for the TCCON FTS retrievals is shown to not add a bias. The same calibration factor is generated using aircraft profiles as a priori and with the TCCON standard a priori. With a calibration to WMO standards, the highly precise TCCON CO2 measurements of total column concentrations provide a suitable database for the calibration and validation of nadir-viewing satellite

Functionalisation of bolaamphiphiles with mononuclear bis(2,2'-bipyridyl)ruthenium(II) complexes for application in self assembled monolayers

A novel ruthenium(II) polypyridyl complex connected competently to a bolaamphiphile, containing amide linkages to provide rigidity via hydrogen bonding in the monolayer, has been prepared. The ruthenium(II) complexes of this ligand and of the intermediates in the synthesis were prepared by modification of the coordinated ligands, demonstrating the synthetic versatility and robustness of this family of complexes. All ruthenium complexes were characterised by electrochemical and spectroscopic techniques and were found to have similar properties to the parent complex [Ru(bipy)[3]][2][+], and remain versatile photosensitisers, with well-defined properties, despite extensive substitution of the bipy ligand

Structural insights into the extracellular recognition of the human serotonin 2B receptor by an antibody

Highly selective monoclonal antibodies recognizing the extracellular 3D epitope of G protein-coupled receptors represent valuable tools for elucidating receptor function and localization in the cell and show promise for a range of therapeutic applications. Here we present the structure of a complex between the human serotonin 2B receptor, captured in an active-like state, and an antibody Fab fragment, bound to the extracellular side of the receptor. The structure uncovers the mechanisms of receptor activation and of extracellular receptor recognition by antibodies

Establishing nonlinearity thresholds with ultraintense X-ray pulses

X-ray techniques have evolved over decades to become highly refined tools for a broad range of investigations. Importantly, these approaches rely on X-ray measurements that depend linearly on the number of incident X-ray photons. The advent of X-ray free electron lasers (XFELs) is opening the ability to reach extremely high photon numbers within ultrashort X-ray pulse durations and is leading to a paradigm shift in our ability to explore nonlinear X-ray signals. However, the enormous increase in X-ray peak power is a double-edged sword with new and exciting methods being developed but at the same time well-established techniques proving unreliable. Consequently, accurate knowledge about the threshold for nonlinear X-ray signals is essential. Herein we report an X-ray spectroscopic study that reveals important details on the thresholds for nonlinear X-ray interactions. By varying both the incident X-ray intensity and photon energy, we establish the regimes at which the simplest nonlinear process, two-photon X-ray absorption (TPA), can be observed. From these measurements we can extract the probability of this process as a function of photon energy and confirm both the nature and sub-femtosecond lifetime of the virtual intermediate electronic state

Crystal structure of rhodopsin bound to arrestin by femtosecond X-ray laser.

G-protein-coupled receptors (GPCRs) signal primarily through G proteins or arrestins. Arrestin binding to GPCRs blocks G protein interaction and redirects signalling to numerous G-protein-independent pathways. Here we report the crystal structure of a constitutively active form of human rhodopsin bound to a pre-activated form of the mouse visual arrestin, determined by serial femtosecond X-ray laser crystallography. Together with extensive biochemical and mutagenesis data, the structure reveals an overall architecture of the rhodopsin-arrestin assembly in which rhodopsin uses distinct structural elements, including transmembrane helix 7 and helix 8, to recruit arrestin. Correspondingly, arrestin adopts the pre-activated conformation, with a ∼20° rotation between the amino and carboxy domains, which opens up a cleft in arrestin to accommodate a short helix formed by the second intracellular loop of rhodopsin. This structure provides a basis for understanding GPCR-mediated arrestin-biased signalling and demonstrates the power of X-ray lasers for advancing the frontiers of structural biology

Anomalous nonlinear X-ray Compton scattering

X-ray scattering is typically used as a weak linear atomic-scale probe of matter. At high intensities, such as produced at free-electron lasers, nonlinearities can become important, and the probe may no longer be considered weak. Here we report the observation of one of the most fundamental nonlinear X-ray–matter interactions: the concerted nonlinear Compton scattering of two identical hard X-ray photons producing a single higher-energy photon. The X-ray intensity reached 4 × 1020 W cm−2, corresponding to an electric field well above the atomic unit of strength and within almost four orders of magnitude of the quantum-electrodynamic critical field. We measure a signal from solid beryllium that scales quadratically in intensity, consistent with simultaneous non-resonant two-photon scattering from nearly-free electrons. The high-energy photons show an anomalously large redshift that is incompatible with a free-electron approximation for the ground-state electron distribution, suggesting an enhanced nonlinearity for scattering at large momentum transfer

Anomalous nonlinear X-ray Compton scattering

X-ray scattering is typically used as a weak linear atomic-scale probe of matter. At high intensities, such as produced at free-electron lasers, nonlinearities can become important, and the probe may no longer be considered weak. Here we report the observation of one of the most fundamental nonlinear X-ray–matter interactions: the concerted nonlinear Compton scattering of two identical hard X-ray photons producing a single higher-energy photon. The X-ray intensity reached 4 × 1020 W cm−2, corresponding to an electric field well above the atomic unit of strength and within almost four orders of magnitude of the quantum-electrodynamic critical field. We measure a signal from solid beryllium that scales quadratically in intensity, consistent with simultaneous non-resonant two-photon scattering from nearly-free electrons. The high-energy photons show an anomalously large redshift that is incompatible with a free-electron approximation for the ground-state electron distribution, suggesting an enhanced nonlinearity for scattering at large momentum transfer